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ABSTRACT 

The Spitzer SDSS GALEX Spectroscopic Survey (SSGSS) provides a new sample of 101 
star-forming galaxies at z < 0.2 with unprecedented multi-wavelength coverage. New mid- to 
far-infrared spectroscopy from the Spitzer Space Telescope is added to a rich suite of previ- 
ous imaging and spectroscopy including ROSAT, GALEX, SDSS, 2MASS, and Spitzer/ 'SWIRE. 
Sample selection ensures an even coverage of the full range of normal galaxy properties, spanning 
two orders of magnitude in stellar mass, colour, and dust attenuation. In this paper we present 
the SSGSS data set, describe the science drivers, and detail the sample selection, observations, 
data reduction, and quality assessment. Also in this paper, we compare the shape of the thermal 
continuum and the degree of silicate absorption of these typical, star-forming galaxies to those 
of starburst galaxies. We investigate the link between star formation rate, infrared luminosity, 
and total PAH luminosity, with a view to calibrating the latter for SED models in photometric 
samples and at high redshift. Lastly, we take advantage of the 5-40 /im spectroscopic and far 



infrared photometric coverage of this sample to perform detailed fitting of the Draine et al. ( 2007 ) 



dust models, and investigate the link between dust mass and star formation history and AGN 
properties. 



Subject headings: infrared: galaxies 
dust, extinction 



galaxies: active — galaxies: ISM — galaxies: stellar content 



1. Introduction 

Spectroscopy in the mid- to far-infrared is 
emerging as a critical tool for understanding the 
physical properties of galaxies, and is the missing 
link in the new era of very large multi-wavelength 
spectroscopic galaxy surveys. This spectral region 
is dominated by emission from interstellar gas and 
dust, largely powered by active star formation. 



Observation of this reprocessed starlight enables 
measurement of embedded star formation that is 
inaccessible to ultraviolet (UV) or optical diag- 
nostics. Spectroscopy in the infrared (IR) also di- 
rectly probes the makeup and physical state of the 
interstellar medium (ISM) through the shape of 
the thermal continuum, aromatic molecular bands 
and a wealth of emission lines. 

The improved understanding of star formation 
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on a galactic scale granted by IR spectroscopy 
is increasingly important as we strive to under- 
stand the evolution of cosmic star formation, and 
in particular its decline since z^l. While IR spec- 
troscopy has enabled substantial progress in map- 
ping the most intensely star-forming (SF) galax- 
ies in the local universe, these galaxies are not 
the principal stellar factories of the current epoch, 
nor are they the closest analogs to the dominant 
sources of star-formation at its cosmic peak. 

At z^l, luminous infrared as well as ultra- 
luminous infrared galaxies (LIRGS and ULIRGS) 
appear to dominate star-forming activity ( |Le| 



Floc'h et al. 20051, in stark contrast to the local 



universe where lower-luminosity 'normal' galax- 
ies contribute the bulk of current star formation. 



However a number of independent studies (Mel- 



bourne, Koo fe Le Flochl[2005l [BelTet al.||2005 



Noeske et al. 20071 Zheng et al. 2007) show that 



these z^l LIRGS and ULIRGS do not display the 
signs of violent interactions observed in their local 
counterparts. In fact, it appears that over much 
of cosmic history the bulk of stars were formed 
in galaxies dynamically similar to normal, disk- 
dominated galaxies at low redshift. 

Thus, low-redshift disk galaxies are important 
laboratories of cosmic star formation — both in 
their own right and as dynamical analogs of higher 
redshift LIRGS and ULIR GS. The Sloa n Digital 
Sky Survey (SDSS) ([York et al.|[2000|) has en- 



abled dramatic strides in our understanding of this 
galaxy population (Bell et al.||2003[ |Brinchmann 



et al. 2004 Kauffmann et al. 



|2003a[ [20041 pair 
ton et al.|2003a||Balogh et al.|2004| ), and GALEX 



observations of the SDSS sample has provided a 
powerful additional lever arm for understanding — 
and selecting for — the star formation properties of 



Yi et al. 



these galaxies (e.g. Heckman et al.[2005 
120051 

Mid-IR (MIR) to far-IR (FIR) spectroscopy 
is now the key missing ingredient. Its addition 
to the GALEX and SDSS data allows a thor- 
ough accounting of star formation, opens up a 
wealth of new diagnostics, and grants the potential 
to self-consistently treat stellar populations and 
dust absorption and emission in models (e.g. Da 
Cunha, Chariot, & Elbaz 2008). A sample of star- 
forming galaxies spanning a representative range 
of physical properties with comprehensive multi- 
wavelength coverage from the far-UV (FUV) to 



the FIR is necessary. 

The Spitzer-SDSS-GALEX Spectroscopy Sur- 
vey (SSGSS) provides such a sample. Selecting 
from the low-extinction Lockman Hole within the 
Sloan+GALEX footprint, FUV and optical diag- 
nostics are used to define a broad, representa- 
tive sample of 100 normal, star-forming galax- 
ies between redshifts and 0.2, boasting rich, 
multi-wavelength coverage. Deep Spitzer low- 
resolution spectroscopy has been obtained for the 
entire sample, and high-resolution spectroscopy 
for the brightest 33 galaxies. 

In this paper we describe the SSGSS dataset. In 
Section [2] we outline the scientific motivation and 
goals of the survey, in Section [3] we describe the 
sample selection and properties, in Section [4] we 
detail the observations, in Section [5] we describe 
the data reduction and quality assessment, in Sec- 
tion [6] we present the spectra and compare to the 
IR spectra of starburst galaxies, and in Section [8] 
we describe the SSGSS data products. 

2. Scientific Motivation 

The primary goal of SSGSS is, broadly, the de- 
tailed characterization of the IR spectra of normal, 
star-forming galaxies spanning a comprehensive 
range of physical properties. Through comparison 
with multi-wavelength data from the FUV to the 
FIR, SSGSS seeks to disentangle the physical pro- 
cesses responsible for the wide range of emission 
features observed in galaxies' IR spectra, and so to 
calibrate these features as diagnostics of galaxies' 
physical states. 

2.1. The Importance of MIR Spectroscopy 

Interstellar dust absorbs from 50% to as much 
as 90% of the UV and optical starlight in star- 
forming galaxies (e.g., Calzetti, Kinney, & Storchi- 
Bergmann 1994; Wang & Heckman 1996; Buat 
et al. 1999, 2002; Sullivan et al. 2000; Bell & 
Kennicutt 2001; Calzetti 2001; Goldader et al. 
2002), and the re-radiated IR photons consitute 
half of the bolometric luminosity in the local uni- 
verse. Broad-band measures of the IR Spectral 
Energy Distribution (SED) of galaxies have taken 
us a long way towards understanding the link be- 
tween this emission and star formation; first with 
the Infrared Astronomical Satellite (IRAS; e.g., 
|Helou||T986l |Rowan-Robinson fc Crawford|[l989l 
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Devereux & Young||1991 


Sauvage & Thuan 


1992 


Buat & Xu||1996| |Walter 


dos & Greenawalt 


199G 


Kennicutt 1998 


Kewley et al. 20021, and 


sub- 



sequently with the Spitzer Space Telescope, with 
deep imaging surveys tracing obscured star forma- 
tion in LIRGs, ULIRGs and normal disk systems 
out to redshift 3 and beyond (e.g. Perez-Gonzalez 
et al. 2005). 

However several processes complicate the inter- 
pretation of broad-band measurements in the IR: 
the thermal IR continuum depends on both the 
temperature distribution and composition of the 
dust; aromatic bands in the MIR are linked to star 
formation, dust composition, and AGN activity, 
but are poorly understood; atomic and molecular 
lines span a wide range of ionizations and hence 
arise from a wide variety of environments; AGN 
and stellar emission compete with thermal dust 
emission in the continuum; and all of these may 
be affected by extinction (although more weakly 
than in optical and UV bands.) To begin to de- 
convolve the effects of these processes, and hence 
to unlock the diagnostic power of the IR spectrum 
of galaxies, high quality spectroscopy is necessary. 

The Spitzer Infrared Spectrograph (IRS) has 
enabled great strides in this regard, especially for 
very luminous IR sources. We now have detailed 
MIR spectroscopic libraries of the most actively 
star-forming galaxies, both in the local universe 



(eg- 


Brandl et al. 2006 Imanishi, Maiolino & Nak- 


agawa 2010 


) , and, increasingly, at higher redshifts 


(eg- 


Yan et al.||2007 Murphy et al.||2009| |Dasyra 


et al.|2009 


Desai et al. 2009; Hernan-Caballero et 


al. 20091. The Spitzer Infrared Nearby Galaxies 



especially important in linking the primary com- 
ponents of the MIR spectrum with physical prop- 
erties, including: linking the effects of metallic- 
ity, ionizing field, star-formation, and H2 inten- 
sity with Polycyclic Aromatic Hydrocarbon (PAH) 



features ( |Smith et aLl|2007| |Roussel et alT1|2007| ) 



physical modeling and characterization of dust 



content (Draine et al. 2007 Munoz-Mateos et al 



2009b), and diagnosing abundances and the rela- 
tive contributions of star formation versus AGN 
activity ( pale et al.||2009||Moustakas et al.||2010 |. 
With its local sample, one of the strengths of 
SINGS is its capacity for spatial resolution within 
its galaxies. This has enabled, for example, the 
mapping of the radial distributions of dust, gas, 



stars, and their evolution ( Munoz-Mateos et al 
2009a|b 2011 ) and the use of Ha attenuation as a 
tracer of obscured star-formation (IPrescott et al 



2007) 



These results have provided a powerful insight 
into the inner workings of galaxies, however it 
is challenging to draw general conclusions on the 
global properties of the SINGS sample due to its 
focus on the nuclear regions of these extended 
sources. In addition, SINGS is typical of other 
IRS samples, in that it is dominated by IR lumi- 
nous galaxies. To characterize the IR spectra of a 
representative population of star-forming galaxies, 
we need increased depth, a greater redshift range, 
and more careful selection criteria. 

2.2. MIR Spectral Diagnostics 

SSGSS spectroscopy spans the 5-40 /xm range, 
covering a wealth of diagnostic features. These 
include: prominent emission bands from PAH 
molecules, a continuum fueled by reprocessed as 
well as direct starlight and by AGN, and abun- 
dant emission lines revealing a very wide range of 
ionization states. Each of these features promises 
significant diagnostic potential. 

2.2.1. PAH Bands 

PAH molecules produce a convoluted spectrum 
of broad features and complexes that dominate the 
MIR below 20 /im. These arise from multiple vi- 
brational modes of aromatic carbon lattices with 
a wide range of grain sizes and ionization fractions 
(see Tielens ( 2005 ) for a review) . 

The SSGSS spectra span the full complement 
of major PAH bands from 6.2 /im to 17 /im, 
and include features arising from carbon-carbon 
stretching, and carbon-hydrogen in-plane and out- 
of-plane bending modes. The strengths of these 
PAH bands, relative to both the continuum and to 
each other, have the potential to serve as powerful 
diagnostics of the environments of these molecules. 
Examples of physical properties probed by PAH 
spectra include: 

Star formation history: the overall luminosity 
radiated in the PAH bands is directly linked to 
the Star Formation Rate (SFR; e.g. Teplitz et 



al. 2007| |Brandl et al.|[2006| ) , and may be an im- 
portant proxy at higher redshifts where PAH fea- 
tures fall in MIR bands. As the PAH spectrum 
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varies depending on star formation history and 
other properties, its diagnostic potential depends 
on a suite of templates measured at low-redshift 
and calibrated to a range of properties more ro- 
bustly measurable at high redshift. 

PAH grain size: grain size distribution affects 
the relative strength of short-to-long wavelenth 
PAH bands. Larger PAH molecules tend to emit 
more efficiently at longer wavelengths than do 



smaller grains ( 


Tielens 


2005 


Schutte et al. 


1993) 


and 



be used as an indicator of grain size distribution. 
This in turn is an indicator of both the growth and 
destruction of PAH molecules, and so will gauge 
the likely governing parameters, such as abun- 
dance, star formation, and AGN activity. 

PAH ionization: ionized PAH molecules radiate 
via carbon-carbon (CC) stretching modes with sig- 
nificantly greater efficiency that do neutral PAH 
molecules, while carbon-hydrogen (CH) modes do 



not show the same dependence (Tielens 2005 1 . 
The ratios of bands arising from CC modes, such 
as those centered at 6.2 /im and 7.7 /im, to CH 
bands, such as the prominant 11.3 /mi feature, 
provide a sensitive indicator of PAH ionization 
fraction, and so probe the FUV photon flux and 
the temperature distribution of the PAH grains. 

As indicated by the models of [Draine &; Li| 
( 2007 ) , a combined suite of PAH ratios can be used 



to disentangle the effects of grain size distribution 
and ionized fraction. 

2.2.2. The Thermal Continuum 

Combined with MIPS 70 /im and 160 /mi flux, 
the SSGSS spectra enable the most comprehen- 
sive characterization of dust temperature distri- 
butions in normal galaxies to date. While the 
60 /mi/100 /im flux ratio has been shown to be 
a good diagnostic of the general shape of the IR 



spectrum (Helou et al. 2000), and hence of dust 



temperature, the MIR continuum can exhibit a 
range of shapes beyond its slope. These varia- 
tions are likely due to variations in the grain size 
and composition of the warm dust component, 
as well as any AGN component. Coupled with 
the UV/Optical dust attenuation curves, spectral 
shapes to 40 /im grant insight into these grain 
properties in galaxies with no significant AGN 
component. Comparison of continuum shapes 



with optical diagnostics of star formation and 
metallicity will be important in disentangling 
these competing effects. 

Furthermore, short wavelength spectra allow a 
better determination of the contribution of ther- 
mal and stellar emission in galaxies over the wave- 
length range from 5 /im to 10 /im, where dust is 
expected to fully dominate the MIR spectrum. 

SSGSS's MIR spectra and FIR photometry al- 
low for a detailed understanding of dust emis- 
sion in these galaxies, and so provide a thor- 
ough account of reprocessed starlight. The multi- 
wavelength nature of this sample is important 
here; optical diagnostics such as the Ha emission 
line and the 4000A break, as well as UV photom- 
etry, give independent measures of star formation 
that may be used to study the efficiency of re- 
processing. This is important for the use of IR 
photometry as a measure of star formation, and 
also to help understand the fate of galaxies' ion- 
izing radiation. GALEX data is especially helpful 
here: the comparison of Lyman-a emission in the 
FUV band to the MIR-FIR intensity and slope 
allows us to investigate the fate of FUV photons, 
and hence escape fraction, as a function of galaxy 
type (see for example Hanish et al.||2010 ). 



2.2.3. Emission Lines 

MIR line ratios provide valuable diagnos- 
tics of HII regions. Electron densities can be 
studied using S III]i 8 .7 /tm/[S IH]33.4 /im and 
[Ne III]i5.5 /tm/[Ne III]36.o /im, while ionization 
ratios can be studied using [Ne III]/[Ne II], [S 
IV] /[S III] and [Ar III]/[Ar II], and elemental 
abundances using various optical and MIR line 
species. 

Many high stellar mass galaxies (log Af* > 10) 
in the local universe and at z=l possess composite 
optical emission-line spectra showing evidence of 
excitation by both young stars and AGN. These 
signatures are degenerate and are necessarily in- 
complete in characterizing reddened AGN. MIR 
ionization indicators such as [Ne III]/[Ne II] are in- 
valuable for decomposing star formation and AGN 
emission because of the minimal levels of attenu- 



ation seen at these wavelengths (Weedman et al 



2006). Comparison of these with other AGN in- 



dicators, such as UV/optical line ratios (Kewley 
et al. 2001 Kauffmann et al. 2003b), and even 
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Fig. 1 — NUV-3.6 fim flux vs. log M* (left) and D„(4000) (right) for the SSGSS sample (solid circles). The 
blue circles are galaxies with 1.5 mJy < F24 < 5 mJy and are observed with only low-res spectroscopy, 
while the red circles have F24 y^m > 1-5 mJy and are observed with both low-res and hi-res modules. Unfilled 
circles are galaxies above the flux and surface brightness cuts but rejected due to high source density in those 
regions of parameter space. Black points are galaxies of the Lockman Hole sample (Johnson et al. 2006), 



the parent sample of SSGSS. Stellar masses are from Brinchmann et al. (2004). 



PAH equivalent widths (EWs), allow us to anal- 
yse the diagnostic power of all of these measures. 
Combined with detailed MIR spectra, we can in- 
vestigate the relative effects of AGN-sourced and 
starburst-sourced hard radiation fields on the ISM. 

The abundance and composition of raw materi- 
als probably has a significant effect on the build-up 
and behaviour of dust grains. To better under- 
stand this effect, MIR abundance indicators com- 
bined with optical measures of gas-phase metallic- 



ity (Tremonti et al. 20041 can be compared with 



PAH bands, thermal dust continuum emission, 
and with MIR silicate absorption. 

The molecular hydrogen IR fluorescence spec- 
trum provides an indication of the physical condi- 
tions in photo-dissociation regions around molecu- 
lar clouds and, while linked to the FUV radiation 
field, may also be produced by shocks in these re- 
gions. 



3. Sample Definition and Properties 

The SSGSS sample consists of 100 star-forming 
galaxies observed with Spitzer IRS in the Lock- 
man Hole region, and is selected to have ancillary 
multi-wavelength imaging and spectroscopy from 
the FUV to the FIR, and to include galaxies with 
SFRs spanning the full range observed in the nor- 
mal, star-forming population, short of the most 
extreme starbursting sources. 

3.1. The Lockman Hole 

The Lockman Hole, situated around a = 10 .7, 
S = 50deg, is a ~10 square degree field with one 
of the lowest column densities of interstellar ma- 
terial in the Milky Way, with N HI < 7 x 10 19 , 
approximately 5% of the Galactic median. This 
makes it ideal for low-surface brightness surveys, 
and so it has been the subject of extensive multi- 
wavelength surveys. These include deep ROSAT 
imaging, deep GALEX (hifuv,nuv 24.5 AB), 
SDSS imaging (m r ~ 22.2 AB) and spectroscopy, 
Deep 2MASS (x6, m K -17.8 AB,|Beichman et al. 
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2003| and Spitzer/SWIRE ([Lonsdale et aDIgOOg]) 



IRAC (m -22.5, 20 for 3.6, 8.0 /im) surveys, as 
well as future programs (e.g. UKIDSS deep sur- 
vey over the full region). The SDSS primary spec- 
troscopic sample targets all galaxies with r < 17.8 
and yields — 100 galaxies deg -1 . 

3.2. Selection Criteria 



The Lockman Hole sample of Johnson et al 



(2006) forms the parent sample for SSGSS, and 
consists of 916 galaxies at z < 0.2 with GALEX- 
SDSS-Spitzer/SWIRE all-band coverage. The 
wide spectral coverage of this sample — from FUV 
to FIR — allows a consistent treatment of dust 
absorption and the corresponding dust emission, 
providing a total accounting of the star formation, 
stellar mass, and dust attenuation for these galax- 
ies. A detailed description of these and other de- 
rived measurements and physical parameters can 
be found in Section pOj SSGSS targets a subset 
of this Lockman Hole sample for observation with 
the Spitzer IRS. These galaxies are selected based 
on two brightness criteria: 

• 5.8 /im surface brightness: Galaxies fill 
the IRS Short-Low slit, and so at short- 
wavelengths we apply a surface brightness 
limit of l5.8f_i m > 0.75 MJy sr _1 , based on 
IRAC 5.8 /im flux and half-light diameter 
(the median/mean half-light diameter is 5 A 
for the entire sample). 

• 24 /im flux: At long wavelengths most of 
the galaxies are unresolved, and so a flux 
limit of F24 Aim > 1.5 mJy is applied. 



The above criteria yield 154 galaxies. We fur- 
ther restrict the sample by performing an even 
sampling of the planes of near-UV (NUV)-3.6 /mi 
vs. stellar mass (from Brinchmann et al. 2004) 
and NUV-3.6 /im vs. D„(4000); in the dense re- 
gions of these spaces, we select the galaxy with 
the highest 24 /im flux. This process ensures 
that the sample spans a comprehensive range of 
physical properties. Figure [l] shows the distri- 
butions of NUV-3.6 /im color versus stellar mass 
and -D n (4000) for the SSGSS galaxies compared 
to the Lockman Hole sample, and illustrates that 
the two span a similar range of galaxy properties. 
The Lockman Hole sample itself is representative 



of the SDSS galaxy sample with GALEX FUV de- 
tections; that is to say, normal, star-forming galax- 
ies, consisting of both early- and late-type, s with 
1.1 < D n (4000) < 2.1. 

The sample observed with Spitzer consists of 
101 galaxies. Due to a fault with the SL spectrum 
for SSGSS 40, the current SSGSS data release con- 
sists of 100 galaxies. We preserve theoriginal in- 
dexing for consistency. 

Within the SSGSS sample, we designate the 
bright sample and the faint sample, with the bright 
sample consisting of the 33 galaxies with F24 Aim > 
5 mJy. The entire sample is the subject of low- 
resolution 5 /im-40 /im IRS spectroscopy, while 
the bright sample is also observed with 10-20 /im 
high-resolution spectroscopy. 

3.3. Sample Properties 

This final SSGSS sample has a redshift range 
of 0.03 < z < 0.2 with median redshift 0.085, and 
a total infrared luminosity (Ltir) range of 3.7 x 
10 9 L Q < L TIR < 3.2 x 10 n i o , with median 3.9 x 
10 10 Lq. The bright sample has redshift range of 
0.04 < z < 0.2 with median redshift 0.077, and an 
Ltir, range of 5.1 x 1O 9 L < Ltir, < 1-9 x 10 10 L© 
with median 6.9 x 10 10 L Q . The flux and surface 
brightness cuts do eliminate galaxies with very low 
stellar masses (2 x 1O 9 M < M < 2 x lO n M ), 
metallicities (8.7 < log(0/H) + 12 < 9.2), and 
extinctions (0.4 < A Ha < 2.3). 

Table [T] details the 101 galaxies in the original 
SSGSS sample, and Figure [2] shows the distribu- 
tion of redshift, stellar mass, Ltir, and D„(4000) 
of the sample. 

[Lit!] 
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Fig. 2. — Distribution of SDSS-measured redshift, Ltir, stellar mass, and D„(4000) for SSGSS galaxies. 
The blue histograms show the 101 galaxies of the entire SSGSS sample, red histograms show the 33 galaxies 
of the SSGSS bright sample, and open histograms show the galaxies of the Lockman Hole sample of |Johnson| 
eTaD pSSel). 
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Table 1 — Continued 



SSGSS 


RA 


DEC 






NUV 




^TIR C L ©) 


M* (M Q ) a 


D„(4000) b 


S/N c 


Ap.Cor. d 


Notes 
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Notes 



1. Bright sample target; observed with SH module 

2. Problematic calibration of SL2 module 

3. Affected by solar flare 

a. Stellar mass from Kauffmann ct al. I 2003a i 

b. £> 7l (4000) from |Kkuflmamt et al.||J0U3a) ' 

c. Signal-to-noise ratio or as^ss low resolution spectra at 22 ^m (sec Scct. |5.3| 

d. IRAC 7.8 fJ.m mag minus synthetic 7.8 fira mag for aperture correction of SL module (see Sect. |5.4| 
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4. New Spitzer Spectroscopy 



4.3. High-Resolution IRS Spectroscopy 



4.1. Instrument, Observing Mode, and 
AORs 

The full SSGSS sample was observed with the 
IRS on board the Spitzer Space Telescope as a 
Cycle-3 Legacy Survey. All spectroscopic observa- 
tions were performed in spectral mapping mode. 
Observations were divided between 29 Astronom- 
ical Observation Requests (AORs), each consist- 
ing of 4 to 8 galaxies observed in cluster mode. 
All galaxies within an AOR were observed con- 
secutively, using guided telescope shifts to move 
between targets. 

4.2. Low-Resolution IRS Spectroscopy 

The entire SSGSS sample was observed with the 
IRS 'low-res' Short-Low (SL) and Long-Low (LL) 
modules. The SL module spans 5.2 to 14.5 /im 
with resolving power R = 60-125 and has slit 
width 3'/6-3"7. The LL module spans 14 to 38 ^m 
with resolving power R = 57 126 and has slit 
width 10'.'5-10 / .'7. Both the SL and LL modules 
are divided into two sub-slits — the first and sec- 
ond orders — abbreviated SL1 (7.4-15.4 /jm), SL1 
(5.2-8.7 /jm), LL1 (19.5-38 fim), and LL2 (14.0- 
21.3 /im). Telescope shifts position the source 
over each order, and so the first and second or- 
ders appears on opposite halves of the detector. A 
third, 'bonus' order from an additional sub-slit is 
obtained with the second order. 

Four offsets in the direction of the slit were used 
to produce a 1-D dither pattern, with two expo- 
sures at each position. This strategy has been 
shown to increase the quality of output spectra, 
in particular for the handling of rogue pixels, un- 
dersampling, and flat-fielding as described in the 
IRS instrument support team document Obtain- 
ing Spectra of Very Faint Sources with the IRS, 



and in Teplitz et al. (2007b 



For low-res spectroscopy our sensitivity target 
was 5cr per resolution element for the continuum 
and PAH band features, after moderate smoothing 
for the faintest sources. To achieve this sensitivity, 
our bright-sample exposure times are 8 minutes in 
both the SL and LL modules and our faint-sample 
exposure times are 8 minutes in the SL and 16 
minutes in the LL module. 



The 33 galaxies in the SSGSS bright sample 
were also observed with the IRS 'hi-res' Short- 
High (SH) module. This module spans 9.9 to 
19.6 /urn with R « 600, and has slit width 4'/ 7. 

For the hi-res spectroscopy our target line flux 
was 1 x 10 -18 W m -2 at 5er, based on the line 
fluxes for [Ne II], [Ne III], [Ar II] and [S III] 
for star-forming galaxies with infrared luminosi- 
ties comparable to those in the bright sample. Our 
hi-res exposure times are 16 minutes, split over 
separate exposures. 

For each hi-res AOR, eight separate sky spectra 
with a total of 16 minute exposure times were also 
taken. 

4.4. Slit Positions 

A crucial aspect of this program is the compar- 
ison of MIR spectral properties to SDSS optical 
properties and so it was essential to ensure over- 
lap between the slit and the 3" SDSS fiber. All 
slits are wider than the SDSS aperture, and so 
the SDSS aperture was centered within each slit 
with at least 0'.'5 accuracy. Acquisition ("peak- 
up") imaging at each AOR enabled us to achieve 
this accuracy, as descrived in Section |4~5| 

Slit orientations were not specified and so are 
random. As the SL and LL observations are taken 
consecutively, the relative orientation of these slits 
on each galaxy is the same as their relative orien- 
tation on the spacecraft (85.9°). SH slit positions 
are randomly oriented with respect to the low-res 
positions. 

4.5. Peak-up Imaging 

To achieve the accuracy needed for our slit po- 
sitioning, we preceeded the observations in each 
AOR with a pointing-mode peak-up image. The 
IRS Peak-Up facility uses centroid positioning 
from an image of the field in the SL detector ar- 
ray to position the science target in the slit. In 
our case, we used bright sources with known off- 
sets to the targets. Despite the high galactic lat- 
itude of our sources, there were sufficient bright, 
low-proper-motion sources in each field to ensure 
quality peak-ups. 

In addition, each galaxy was observed with ded- 
icated 4x30 second peak-up imaging. These peak- 
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Fig. 3. — Example SL1 (upper) and SL2 (lower) frames through the reduction process. Left is the post- 
pipeline combined frame, middle is the sky-subtracted frame, and right is the cleaned sky-subtracted frame, 
from which 1-D spectra are extracted. Pointers mark the ends of the spectra — to the left of the frames for 
SL1, and to the lower right for SL2. The segments at the upper right are the bonus order spectra. 



ups utilized the blue filter, with spectral coverage 
of 13.3-18.7 /im. This ~ 16 /im coverage gives us 
photometric data intermediate to the IRAC 8 /jm 
and MIPS 24 /jm imaging. 



Figure 18 (Appendix A) shows the dedicated 
peak-up images for all galaxies in the sample, in- 
cluding slit positioning. 

4.6. Solar Flare 

A solar flare during one AOR resulted in very 
high cosmic ray counts for five targets. The af- 
fected galaxies are SSGSS 20, 22, 23, 24, and 25. 



5. Data Reduction 

5.1. Low-Resolution Spectroscopy 

Standard IRS processing was performed by the 
Spitzer Pipeline version S15.3.0. Pipeline steps 
included ramp fitting, dark subtraction, droop, 
linearity correction and distortion correction, flat 
fielding, masking and interpolation, and wave- 
length calibration. 

No separate sky images were taken with the 
low-res modules and so sky frames were con- 
structed from combined science frames. Between 
the first- and second-order exposures, the source is 
shifted to opposite sides of the detector. This re- 
sults in a sky exposure for the order not currently 
containing the source. Within each AOR, all such 
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5 10 15 20 25 

Rest Wovelength (/im) 



Fig. 4. — Example of method used to stitch the low-res spectra. The colored dashed spectra are from 
individual channels, from left to right: SL2 (blue), SL1 (red), LL2 (blue), LL1 (red). The dark gray spectrum 
is stitched from these channels. Dark gray shaded areas represent the channel-overlap regions where two 
channels are medium-combined. The light gray shaded areas are the channel-overlap regions where the 
stitched spectrum is set to the values of the nearest channel. The default medium-combined region is the 
central 50% of the overlap region, however this is adjusted manually to avoid end-of-channel transmission 
drops. The dotted line is the stitch point for the SL1-LL2 junction, which defaults to the middle of the 
overlap region but in many cases is set manually due to transmission drops. 



sky exposures were median combined to produce 
a master sky frame for each order. 

Bad and rogue pixels masks for each exposure 
were created with IRSCLEAN (version 1.9), with 
AGGRESSIVE keyword set to 2, and with nega- 
tive pixels also masked. We also applied the cam- 
paign rogue pixel masks released by Spitzer. 

Figure [3] shows example frames for the SL mod- 
ule through the reduction process, including the 
initial post-pipeline frame with exposures com- 
bined, the sky-subtracted frame, and the cleaned 
frame. 

Extraction of 1-D spectra was performed us- 
ing SPICE version 2.0.1. The default setting 
was used to create a wavelength collapsed spa- 
tial profile, and the center of the profile was de- 
termined manually due to excess failures of the 
auto-center option. SPICE then performed spec- 
tral extraction along this peak. The point source 
tune module applied flux corrections for changes 
in the wavelength-dependent Point Spread Func- 



tion (PSF) width and aperture size, and assumed 
emission by a point source. 

5.1.1. Low-res Combining and Stitching 

The extracted 1-D spectra were combined by 
weighted mean. First, the spectra at the four 
dither positions along the slit were were combined 
for each order. Then the first, second and bonus 
orders were stitched to produce separate SL and 
LL spectra using the following semi-manual pro- 
cess: by default, the central 50% of each overlap 
region was combined by weighted mean and the 
outer quartiles of the overlap region were set to 
the nearest order. In select cases the combined re- 
gion was adjusted to avoid obvious drops in trans- 
mission at the ends of either order. The SL and 
LL modules were combined by a clean cut between 
the spectra, determined manually to avoid end-of- 
channel transmission drops. Figure [4] illustrates 
the method used for stitching spectra. 

The overlap region of the SL and LL mod- 
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ules is ~ 14-14.5 /tm, and in this region the slit 
widths are 3'.'7 and 1C//5 for SL and LL respec- 
tively. The approximate LL-to-SL ratio of admit- 
ted light from centrally-aligned point source in the 
overlap region is ~1.3, based on a simple model 



of the Spitzer PSF (Smith 2004), and the over- 



all LL sensitivity is a factor of ^5 higher than the 
SL module in this region (Spitzer Observer's Man- 
ual, p. 160). While these differences are accounted 
for by the pipeline reduction assuming a centrally- 
aligned point source, for extended sources we may 
expect an offset at their intersection when the 
spectra are stitched. 

In fact no galaxies exhibit noticeable disconti- 
nuities across the modules in this region, indicat- 
ing that deviations from the point-source approx- 
imation in the overlap region are small, even for 
the SL module at ~14-15/«ti. No re-scaling was 
necessary or performed. 

5.2. High-Resolution Spectroscopy 

As with the low-resolution spectroscopy (Sect. 
5.1), the Spitzer Pipeline version S15.3.0 was used 
to perform standard IRS processing. Unlike the 
low-res data, for the hi-res observations we com- 
bined the eight separate exposures of each galaxy 
before extraction. This improved the signal to al- 
low higher quality extractions. 

As noted in Section |4.3[ separate hi-res sky 
frames were taken for each AOR, median com- 
bined from eight 16 minute exposures. Masking 
followed the same method used for the low-res 
spectroscopy, utilizing IRSCLEAN (version 1.9) 
and the Spitzer campaign rogue pixel mask. Like- 
wise, extraction was performed as for the low-res 
data, using SPICE v. 2. 0.1. The hi-res orders were 
stitched into a single spectrum using a similar 
method to the low-res spectra: order-overlap re- 
gions were combined by taking the weighted mean 
within the central 50% or the overlap regions, 
while for the outer quartiles the nearest order was 
used. 

Figure [5] shows example frames for the SH mod- 
ule through the reduction process, including the 
initial post-pipeline frame with exposures com- 
bined, the sky-subtracted frame, and the cleaned 
frame. 
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Fig. 6. — Signal-to-noise ratios of low-res spectra 
at 10 fim (upper) and 22 /im (middle) versus IRAC 
7.8 fan and MIPS 24 fan mag respectively, and of 
the hi-res spectra at 15 fan versus MIPS 24 /mi 
mag. Dashed lines indicate target S/N. 



5.3. Achieved Sensitivity 

We estimate the continuum sensitivity achieved 
by determining the signal-to- noise ratio (S/N) in 
the processed spectra at relatively featureless re- 
gions of the spectra. We compare the mean flux 
to the corresponding noise in the regions 9-10 /mi 
and 20-24 fan for the low- res spectra and 14.7- 
15.3 /im for the hi-res spectra. The noise is taken 
to be the standard deviation in these regions af- 
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Fig. 5. — Example SH frames through the reduction process. Left is the post-pipeline combined frame, 
middle is the sky-subtracted frame, and right is the cleaned sky-subtracted frame, from which 1-D spectra 
are extracted. 
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Fig. 7. — Signal-to-noise ratios of the [Ne II] (up- 
per) and [Ne III] (lower) emission lines from the hi- 
res spectra versus the PAHFIT-measured fluxes. 
Dashed lines indicate target S/N. 



ter subraction of a power law fit. Figure [6] shows 
the continuum S/N for low-res and hi-res spec- 



tra versus IRAC 7.8 fim and MIPS 24 /im mag. 
At ~ 22 /im nearly all sources exceed the target 
S/N = 5, with a minimum S/N of 4, and have 
the expected tight relationship with 24 /im mag. 
At ~ 10 /im the variable silicate absorption re- 
sults in a number of sources with S/N < 5. In this 
case there is greater variability of S /N with photo- 
metric magnitude due to this absorption and the 
overlap of more PAH bands with the IRAC 7.8 /im 
band. 22 (im S/N is provided in Table [l] 

We estimate the S/N for the [Ne II] 12.8 and 
[Ne III] 15 .6 emission lines in the hi-res spectra by 
comparing the integrated flux of the lines to the 
standard deviation in the relatively smooth, blue- 
ward regions at 12.25-12.75 /im and 15-15.5 /xm. 
Figu re [7 shows line S/N versus the line flux (see 



Sect. 6.31. [Ne II] and [Ne III] are detected in all 



galaxies with S/N > 5. In most cases, both [Ne II] 
and [Ne III] are detected with S/N > 10. There is 
strong scatter in the trend between line flux and 
S/N due to the weakness of the trend between con- 
tinuum and line luminosity, and this is especially 
pronounced in [Ne III], where the line emission 
may be linked more to an AGN components than 
to star formation. 

The flux uncertainties reported by SPICE, 
which are included in the released data, have been 
checked for the low-res modules by comparing the 
multiple dithered observations for each galaxy. 
These uncertainties were found to be accurate to 
within a factor of 1.5. 
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5.4. Aperture Effects 




Fig. 8. — Synthetic minus SWIRE photometry as 
a function of r-band petrosian radius 8 /im, 16 /im, 
24 /im. The trend between this magnitude differ- 
ence and radius observed in the 8 /im band sug- 
gests slit losses in the SL module, but not at longer 
wavelengths. 

An increasing fraction of the wavelength- 
dependent Spitzer PSF falls outside a given IRS 
slit with increasing wavelength. The standard 
spectral extraction with SPICE corrects for this 
slit loss assuming emission by a point source. At 
> 16 /im the combination of the larger Spitzer PSF 
and the wider IRS slits means that all galaxies 



are well-approximated as point sources (Fig. 18 
Appendix A). However, a number of the SSGSS 
galaxies have sufficient spatial extent that they 
are not well-approximated by point sources in the 
narrower SL and SH modules. 

We investigate slit loss in the low resolution 
modules by comparing synthetic photometry from 
the reduced, calibrated low-res spectra to large- 
aperture photometry obtained from SWIRE imag- 



Fig. 9. — SL slit aperture magnitude minus syn- 
thesized and circular aperture magnitudes as a 
function of r band petrosian radius. 



ing (IRAC channel 4 7.8 /im and MIPS 24 /im), 
and from the SSGSS 16 /im peak-up imaging. Syn- 
thetic photometry is obtained by integrating the 
spectra using k_correct (v4.14; Blanton03b). Fig- 
ure [8] shows the difference in these magnitudes as 
a function of the r' band Petrosian radius [Rp), 
with aperture photometry performed in 7" and 
12" radius apertures for 8 /im and 24 /im respec- 



tively, as described in Johnson et al. (2007). Dif- 



ferences may be attributed to a number of sources 
(e.g., source confusion in the imaging aperture 
photometry, relative calibration errors) but any 
trend with the angular size of the galaxy is sug- 
gestive of slit losses. No such trend is observed at 
16 /tm and 24 /im. However, at 8 /im a significant 
trend is apparent, indicating slit losses in the SL 
module. 

To isolate the effect of slit loss from any calibra- 
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tion uncertainties, we perform photometry of the 
8 /im IRAC images within the region correspond- 
ing to the SL slit aperture (a 3''65x9?6 rectan- 
gle), and compare this to the synthetic and circu- 
lar aperture magnitudes. Figure [9] (upper) shows 
the difference between slit aperture and synthetic 
magnitudes. A 0.3 mag aperture correction is ap- 
plied by SPICE at 8 /im and so the synthetic mag- 
nitudes are higher than the derived slit aperture 
magnitudes. Taking this offset into account, the 
slit aperture samples more flux than the extracted 
spectra, with an average offset of ~ 0.2 mag and 
similar scatter. The magnitude difference is inde- 
pendent of galaxy size, indicating that the trend 
observed in Figure [Sldoes not result from calibra- 
tion effects. Figure|9] (lower) shows slit aperture 
versus circular aperture magnitudes, revealing a 
similar trend to that observed with the 8 /im syn- 
thetic magnitudes. The maximum difference is 1- 
1.5 mag for R p {r') > 9", 

The apparent slit losses in the SL module mean 
that care is necessary when considering the shape 
of the (< 16 /im) Spectra Energy Distribution 
(SED) of large galaxies. The difference between 
the synthetic and photometric magnitudes at 8 /im 
should provide a reasonable estimate of the aper- 
ture correction. These corrections are provided 
Table [U 

6. SSGSS Data 

Figure [18] shows multi- wavelength imaging for 
all 101 galaxies in the SSGSS sample, including 
GALEX FUV and NUV, SDSS multi-colour, r, 
H, IRAC 3.6 /im and 8 /im, 16 /im from SSGSS 
peakup images, and MIPS 25 /im and 70 /im. 
Overlaid are the slit positions for the SL, LL, and 
SH modules. 

6.1. Low-Resolution Spectra 



Figure 19 (Appendix B) shows the reduced, 
combined, and stitched SL and LL spectra for 
the 100 galaxies in the current data release. For 
each galaxy the SSGSS catalog number, BPT type 
(paldwin et al. ||1981| |Kewley et ai~1|2001| [KailfF 
mann et al.|2003b ), redshift, stellar mass, NUV-R 
color, and -D„(4000) are given. 

The relatively low sensitivity of the 5.25- 
7.58 /im SL2 module, coupled with the relative 
faintness of many galaxies in this range, results 



in very low SL2 signal in 10 sources. This spec- 
tral region should be used with caution in these 
sources, which are noted in Table [TJ 

For one galaxy, SSGSS 56, the SL-LL over- 
lap region falls directly on the redshifted 12/im 
PAH complex and the [Ne II] 12.8 line. End-of- 
order transmission drops in both the SL and LL 
modules affect these features, and so both are lost 
in this spectrum. 

6.2. High-Resolution Spectra 



Figure 20 (Appendix B) shows the reduced, 
combined, and stitched SH spectra for the 33 
galaxies in the SSGSS bright sample. For each 
galaxy the SSGSS catalog number, BPT type, red- 
shift, stellar mass, NUV-R color, and L>„(4000) 
are given. 

6.3. Emission Line Strengths 

We measure emission line strengths by integrat- 
ing line regions after subtraction of the continuum. 
Here, the continuum refers to all non-emission line 
features, including PAH features and the thermal 
continuum. The continuum is determined using 
the PAHFIT code of iSmith et all d2007b, which 



fits a physically-motivated model to each spec- 
trum using x 2 minimization. We describe our 



use of PAHFIT in detail in |Q'Dowd et al.| (|2009j) 
(OD09). When applied to the hi-res spectra, we 
remove emission line fitting from the process and 
mask the emission line regions. We are not con- 
cerned with accurate measurement of continuum 
features — only that we find a reasonable extrapo- 
lation of the continuum in the region of the emis- 
sion lines. 

The emission line strengths measured by PAH- 
FIT from the hi-res spectroscopy are highly reli- 
able, as the continuum near most lines is reason- 
ably smooth, allowing robust subtraction. These 
line strengths compare closely to those measured 
using a simple linear fit to the continuum sur- 
rounding each line. In the case of the [Ne II] 12.8 
line, PAHFIT is significantly more reliable than 
cruder continuum-fitting methods because of the 
close proximity of this line to the peak of the 
12.7 /im PAH feature. 

Figure 10 compares [Ne II] 12.8 an d [Ne III] 15.6 



lines strengths between the low- and hi-res spec- 
tra. Although there is general agreement, on av- 



16 



1CT 



E 
5 



g 1 10" 



in 



c 



10" 

o 

u 
0) 



4> 



10" ,s - 




10"'° 10" 
Low-res Integroted Line Strength (Wm -2 ) 



10~' 6 



Fig. 10. — Comparison of the [Ne II] 12.8 A™ (blue) and [Nc III] 15.6 /im (red) emission line integrated strengths 
as measured with PAHFIT using the low-res and hi-res spectra. Error bars represent the 1-a uncertainties 
in the fits. 



erage the low-res spectroscopy shows weaker line 
strengths by ~20-50%, with similar scatter com- 
pared to the hi-res spectra. In four galaxies, PAH- 
FIT claims detection of [Ne III] lines in the low-res 
spectra where they are undetected at high resolu- 
tion. The uncertainties reported by PAHFIT for 
the low-res [Ne II] lines are significantly underes- 
timated, suggesting that there is more degeneracy 
between this line and the 12.6 /imPAH complex 
than estimated by in the fit. 

Table [2] gives emission line strengths measured 
from the hi-res spectra. Line strengths from the 
low-res spectra can be found on the SSGSS website 
(see Sect. 



• 2) 



6.4. Composite Spectra 



Figure 



11 



(upper) shows the composite of the 
low-res SL and LL spectra for the entire SSGSS 
sample, normalized between 9-11 /im. Primary 
PAH features are prominent, including a pro- 
nounced 17 /im complex. Several emission lines 
are also distinct. 

(lower) shows the composite of the 



Figure 
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SH spectra for the SSGSS bright sample, normal- 
ized between 14-15 /an. PAH features have sig- 
nificantly more detailed line profiles than in the 



low-res spectra, including clear resolution of the 
16.5 /im feature from the 17 /im complex and inter- 
esting structure in both the 11.3 /im and 12.7 /im 
features. In addition, the [Ne III]i5.6/im line is 
well-resolved from the 12.7 /im feature. 

7. Results and Discussion 

7.1. Comparison to Starburst Spectra 

The SSGSS sample spans the range of physical 
properties of 'normal' star-forming galaxies, as de- 
scribed in Section [3| We investigate the difference 
between the composite SSGSS spectra and galaxy 
spectra of previous IR surveys which have typi- 
cally been dominated by more extreme starburst 
galaxies. 

Figure [12] shows the low-res SSGSS composite 
spectrum (normalized between 9-11 /im for com- 
posite construction) overplotted with an ISO spec- 
trum of M82 dSturm et al||2000h and the Spitzer 



local starburst composite of Brandl et al. ( 2006 1 



These spectra are normalized according to the re- 
sults of the fits described below. Also plotted is 
the 0.8 < z < 1.2 LIRG composite of |Dasyra et| 
al. ( 20091) , normalized at the blue end of this spec- 



trum, 6.5-6.8 /im. 
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Table 2 
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Fig. 11. — Upper panel: composite SL and LL low-res spectrum of the 100 galaxies of the current SSGSS data 
release (solid black line). Low- res spectra are normalized between 20-24 fxm. The shaded region represents 
1-sigma scatter of the normalized spectra about the mean spectrum. Lower panel: composite SH hi-res 
spectrum of the 33 galaxies of the SSGSS bright sample. Hi-res spectra are normalized between 14-15 /xm. 



M82 and the local starburst composite have 
similar EWs for PAH features at 11.3 /im and blue- 
ward. However the EW of the 12.7 /im feature is 
lower in M82, and the EW of the 17 /zm feature is 
lower in both. These starburst galaxies appear to 
have stronger long wavelength continuum emission 
than the galaxies of the SSGSS sample. Note that 



the MIR spectral shape of M82 is within the range 
seen for the Brandl et al. (2006) starbursts. The 



composite for z ~ 2 ULIRGs shows significantly 
reduced EWs for all PAH bands, and a thermal 
dust component slope consistent with that of the 
local starburst template, however the short wave- 
length range of the ULIRG composite makes this 
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Fig. 12. — Top: Composite low-res spectrum of the SSGSS sample (solid black line), with the spectrum of 
M82 (green; Sturm et al 2000), the composite spectrum of local starbursts (blue; Branch et al. 20061, and 



the composite of z ~ 1 ULIRGS (magenta; Dasyra et al. 2009). The shaded region represents the 1-sigma 



scatter of the normalized spectra about the mean spectrum. The dotted black lines show best fits to M82 
and starburst composites of a model consisting of the SSGSS composite plus a dual-temperature blackbody. 
The fits indicate that starburst spectra are well represented by that of a normal, star-forming galaxy plus 
an enhanced thermal dust component at T < \2QK. Bottom: 10-20/im region of top panel. 



difficult to establish from these data. 

To determine whether the difference in spec- 
tral shape between local starbursts and SSGSS 
galaxies is consistent with a difference in small- 
grain dust temperature distribution, we fit the 
M82 and local starburst composite spectra with 
a model consisting of the SSGSS low-res compos- 



ite spectrum plus a dual-temperature blackbody 
dust component, in the case of M82 masking a 
10 /im window around the 9.7 /^m silicate absorp- 
tion feature. Both are well fit with the addition of 
a 120K component plus a single cooler component 
at around 80K, and are poorly fit with components 
warmer than 130K. The lower EWs of the 12.7 /jm 
and 17 /xm features are reproduced in the fit to the 
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starburst composite, while only the 12.7 /im fea- 
ture is accurately reproduced in M82; the 17 /zm 
feature has lower EW than can be explained by a 
heightened continuum. In OD09 it was seen that 
the 17 /im feature is relatively weak compared to 
the 7.7 /im feature in galaxies undergoing more in- 
tense star formation. This is consistent with our 
observation of the weaker feature in M82. The re- 
sults of these fits are also shown in Figure [12] The 
relative normalization of all spectra in this figure 
are determined by the best fits. 

If we assume that the emission between 20 /im 
and 30 /jm in the SSGSS composite arises solely 
from a 120K thermal dust continuum component, 
then we estimate the additional flux from this 
component as a factor of two higher in the local 
starburst composite and a factor of eight higher 
in M82 compared to the SSGSS average. Al- 
though the additional continuum light apparent 
at the longer wavelengths of the starburst spectra 
is likely to be due to dust at a range of temper- 
atures, this simple model shows that the differ- 
ence in spectra is consistent with starburst galax- 
ies having a higher proportion of emission from 
< 120K dust than more quiescent galaxies. 

7.2. Comparison of Absorption and Con- 
tinuum Properties 

The shape of the MIR spectra of the most ex- 
treme IR emitters tends to be dominated by either 
a bright continuum — from either AGN or hot 
dust — or by heavy silicate absorption centered 
at - 9.7 (Mm. |Spoon et ai] ( |2007[ ) (S07) showed 
that this gives rise to a bimodality in the space 
of 6.2 /jm EW versus silicate absorption strength 
for ULIRGs and other IR-luminous galaxies and 
for AGN, from analysis of sources observed with 
Spitzer and ISO. 

Figure [13] shows the location of the SSGSS 
galaxies in this space, with the regions occupied by 
the S07 sources shaded. We calculate our indica- 
tor of absorption strength similarly to the method 
used by S07 for starburst galaxies. We take the 
ratio of the 9.5-10.5 /jm minimum to the "unab- 
sorbed" flux at this point, the latter determined by 
interpolation using a linear fit between 4.5-5 /im 
and 14.5-15 /zm regions, which are relatively un- 
affected by PAH emission or silicate absorption. 
6.2 /im PAH EW is calculated using PAHFIT 
to measure the feature strength, as described in 



Treyer et al] ( |2010[ ) (hereafter T10). This method 
differs from that used by S07, who determine the 
continuum, and hence the EW, by fitting a spline 
between the edges of the 6.2 /im feature. 

Methods for EW determination based on esti- 
mation of the continuum from the boundaries of 
the 6.2 /im feature result in saturation at widths 
greater than ~0.5-l /im. This is due to the fact 
that, at higher line strengths, the continuum esti- 
mation includes increasing amounts of light from 
the 6.2 /zm feature and the neighboring 7.7 /im 
PAH complex in the continuum estimation. At 
EW > 0.5-1 /im the flux at the points typically 
used to estimate the continuum is in fact domi- 
nated by PAH emission. Sargsyan fc Weedman| 
([2009]) report 6.2 /zm EWs for SSGSS galaxies, 



finding EW < 1 /im for the entire sample. They 
calculate EW using a continuum fit anchored at 
5.5 /im and 6.9 /im. This approach, like the spline 
method, is vulnerable to underestimates and sat- 
uration at high EWs due to the presence of sig- 
nificant PAH flux in the regions used for contin- 
uum estimates. T10 presents an analysis of this ef- 
fect, and a comparison with the [Sargsya n fc Weed-| 
(2009) results. As many SSGSS galaxies have 



man 



strong 6.2 /im features, it is important that their 
EWs be measured using full continuum+PAH pro- 
file fitting. 

In Figure [13] we can see that SSGSS galaxies 
have silicate absorption indicators consistent with 
the low absorption locus. A number of our sources 
are consistent with the starburst galaxies in S07's 
analysis, but much of our sample extends to sig- 
nificantly higher EWs, in part because of the EW 
saturation issue. Sources with AGN components 



by BPT designation ( 


Baldwin et al. 1981 


Kewley 


et al.|2001 Kauffmann et al.|2003b 


) have brighter 



continua, and hence lower EWs, than purely Star- 
Forming (SF) galaxies. These EWs may also have 
been reduce by destruction of PAH by the AGN. 

7.3. PAH Contribution to the IR Lumi- 
nosity 

For star-forming galaxies, total infrared lumi- 
nosity {Ltib) is an good proxy for the rate of 
current star formation, with much of the starlight 
in dusty starburst environments being reprocessed 
into the IR. SSGSS has already shown that PAH 
emission is also strongly tied to a galaxy's star 
formation history; the relative intensities of PAH 
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Fig. 13. — Equivalent width of the 6.2 /im PAH feature versus a silicate absorption strength, defined as 
the log of the ratio of the minimum in the 9.7 /im absorption trough to the interpolated value at this point 
from a linear fit to the 4.5-5 /im and 14.5-15 /im regions. SSGSS galaxies are colour-coded according to 
BPT designation ( Kewley et al.pOOT Kauffmann et al. 2003b); blue points are star- forming galaxies, green 
points are composite sources, and red points are AGN. The shaded region shows the positions on this plot 
of LIRGs, ULIRGS, starbursts and AGN from Spoon et al. (2007). These more extreme galaxies tend to 



have stronger silicate absorption and brighter continua than the more quiescent star-forming galaxies of the 
SSGSS sample. 



bands trace specific star formation indicators 
(OD09), while the luminosities of individual PAH 
bands correlate strongly with total SFR (T10). 

Understanding this latter relation, and in gen- 
eral the relative strength of the PAH contribu- 
tion to galaxies' MIR emission, is important for 
modeling the IR SEDs of higher redshift galax- 
ies where the PAH bands cannot be so clearly re- 
solved, and for SED modeling in large photomet- 
ric surveys such as SWIRE ( |Lonsdale et al.|2003| ). 
T10 demonstrated the near-linearity of the rela- 
tion between the intensity of individual PAH fea- 
tures and both Ltir and SFR. For the purpose 
of model fitting to MIR data in which individual 
PAH features are not well resolved, it is also valu- 
able to determine the link between total PAH lu- 
minosity (Lpah) and SFR and Ltir- 



We measure Lpah using PAHFIT (|Smith et 
see Sect 



al. 2007 



6.3), summing the integrated 



strengths of all PAH bands between, and includ- 
ing, the 6.2 /mi to 17 /im features. The rela- 
tively weak 3.6 /im feature falls outside our spec- 
tral range and is not included in our determination 
of Lpah- We correct Lpah for aperture effects, 
comparing SWIRE and SSGSS peak-up photome- 



try to synthetic photometry (see Sect. 5.4 ), and in- 
terpolating corrections to PAH peak wavelengths. 



Ltir is determined by fitting Draine & Li (2007) 



model SEDs to the Spitzer (IRAC, IRS blue peak 
up, and MIPS) photometry, and then integrat- 



ing between 3 and 1100 /im. SFR is from |Brinch- 
mann et al. (2004), who fit the models of IChar- 



lot & Longhetti (2001 1 to strong emission lines of 



galaxies in the SDSS spectroscopic sample, using 
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Fig. 14. — Left: total PAH luminosity versus SFR e . Middle: as with the right panel, but, only galaxies 
with SSFR < 1 x 10~ 10 yr~* are filled, revealing an especially tight trend for galaxies dominated by young 
populations. Dashed lines show the best regression fits to star-forming galaxies (left) and the filled circles 
(middle). Right: standard deviation in the fractional residuals in infrared luminosity after subtraction of 
the best linear fit to the plot of Lpah SFR e . For each point in this plot, the fit was made to the subset 
of the sample with specific star formation indicator (SSFR, D„(4000), or Ha EW) below the given value. 
The trend indicates that, while young populations exhibit a very tight relationship between PAH or total 
IR luminosity and SFR, this relationship becomes increasingly scattered as populations age. 



a Kroupa IMF (Kroupa 2001) and the Chariot 



& Fall (20001 dust model. We use the aperture- 



corrected median of the derived SFR likelihood 
distributions, designated SFR e . 



Figure 14 (left) shows Lpah versus SFR e . 
These quantities are strongly correlated for the full 
SSGSS sample, however we are interested in cal- 
ibrating Lpah to star formation, and so we per- 
form regression analysis for SF galaxies only: 



logw Lpah(L q ) — 
[0.99 ± 0.05] log w SFR e (M Q yr^ 1 ) + [8.8 ± 0.1] (1) 

The slope consistent with that presented in T10 
for the 7.7 ^m complex. As also found by T10, 
galaxies with AGN components generally fall on 
the same locus as purely star-forming galaxies, al- 
though have significantly greater scatter. This is 
primarily due to the fact that SFR is determined 
with much lower accuracy for the more evolved 
populations in which BPT-designated AGN are 
exclusively found. 

This is further explored in Figure [14] (middle) , 
which highlights the youngest 50% of the SSGSS 



sample, as defined by specific SFR e (SSFR; SSFR 
> 10~ 10 yr _1 ), chosen simply to divide the sam- 
ple in two. The correlation is significantly tighter 
when only these more actively star-forming galax- 
ies are included. Fig. 14 (right) shows the stan- 



dard deviation in the log residuals for Lpah af- 
ter subtraction of the best linear fit for a range 
of subsample sizes. The subsamples are sorted to 
contain the youngest stellar populations, selected 
by three different measures: SSFR, D„(4000), and 
Ha EW. In all cases, scatter decreases as the dom- 
inance of the young population increases; almost 
monotonically for D„(4000) and Ha EW, and for 
these measures the scatter flattens significantly for 
the subsample containing the youngest ~60% of 
sources, and for smaller "youngest" subsamples. 
The diagnostics corresponding to the ~ 50 th - 
60 th percentile transition are L>„(4000) = 1.4, 
Ha EW = 18A, and SSFR = 10~ 10 



yr 



For 



more evolved sources, the exceptionally tight trend 
between log Lpah and log SFR e (a < 0.2) gains 
significant scatter. Again, this is primarily due to 
the inaccuracy with which SFR is measured for 
evolved populations. 

(upper panels) show the mass- 



Figure 15 
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Fig. 15.— Upper: L PAH /M„ versus D„(4000) (left), Ha EW (middle), and SSFR (right). Galaxies are 
colour-coded according to BPT type as in Figure [T3| AGNs and composite sources sit on the same loci as 
SF galaxies, and have systematically older stellar populations. Dashed lines show regression fits for the SF 
population (blue points), with AGNs/composites excluded. Lower: As above, but for Lpah / 'Ltir- Dashed 
and dotted lines show the median values for SF galaxies and AGNs/composites respectively. 



weighted relationships. The ratio of Lpah to 
stellar mass (M* is plotted versus optical star for- 
mation diagnostics from SDSS: th e narrow band 



4000A break measure D„(4000) (Kauffmann et 
aLl|2003a[ ) and Ha EW (left and middle respec 



tively), and also SSFR (right), which is described 
above. There are tight linear correlations be- 
tween all quantities, although again we see the 
increased scatter in SFR for evolved galaxies. The 
best-fit regression lines for the SF galaxies are 
also shown. The mean 1-a residual scatter in 
Lpah/M* for SF galaxies after subtraction of 
the best-fit lines is ^20% for all three cases. For 
L>„(4000) and Ha EW, the residual is -20% even 
including galaxies with AGN components, indicat- 
ing that these age diagnostics are good predictors 
of mass- weighted PAH intensity for most galaxies. 



This emphasizes that the increased scatter seen 
in Lpah and Ltir versus SFR in evolved pop- 
ulations is due to inaccuracy in measuring SFR 
rather than, say, a blurring due to evolved stars 
contributing more significantly to dust heating. 



As noted by Calzetti et al. (2010), the correlation 



between PAH emission and metallicity leads to 
another source of increased uncertainty in the link 
between star formation and IR dust emission in 
evolved sources. 

Total IR luminosity is also closely related to 
these SF diagnostics, and has a very similar slope 
and scatter as Lpah, for ah three parameters and 
for both the absolute and mass- weighted relations. 
For the purpose of modeling high-z and photo- 
metric SEDS, it is useful to understand how the 
PAH-to-total IR ratio changes with these proper- 
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ties. Figure [l~5| (lower panels) show Lpah/Ltir 
versus these same SF diagnostics, with dashed and 
dotted lines showing the medians for the SF and 
AGN/composite subsamples respectively. In all 
SF diagnostics there is a significant but highly 
scattered trend, with more evolved galaxies and 
AGN showing weaker PAH features for a given 
L TIR . It is not possible to extricate the effects of 
age and AGN status on L PAH /L TIR from these 
data, although as shown in OD09 both seem to 
play a part. More specifically, short wavelength 
PAH features (6.2 /im and 7.7 fim) are dimin- 
ished in more evolved galaxies and AGN, while 
longer wavelength features are not. This is seen 
in the PAH ratios (OD09) and with respect to 
Ltir (T10). The trend is between Lpah/Ltir 
and star formation is still weakly significant for 
SF galaxies alone, with AGN excluded, although 
this is entirely driven by a small number of heavily 
star- forming sources with very high Lpah/Ltir- 

7.4. Dust Mass 

The SSGSS IR spectra, with its coverage of all 
major PAH features, combined with MIPS 70 /jm 
and 160 ^m photometry, enables us to fit detailed 
dust models to determine dust mass. We use the 
models of Draine & Li (2007), which calculate 



ilar to the D10 SDSS sample. The SINGS dust 
masses are roughly consistent with the range of 



IR emission spectra for starlight-heated dust con- 
sisting of amorphous silicate and graphite grains, 
and of PAH grains. SSGSS spectra and pho- 
tometry are fit to template libraries derived from 
these models (Draine & Li, private communica- 
tion) to estimate the distribution of illuminating 
starlight intensity and the relative abundance of 
PAH. These are combined with MIPS fluxes and 



the relations detailed in Draine & Li ( 2007 ) to de- 



termine the total dust mass. 

Based on these models, the median dust mass 
for the sample is M dust = 8.5 x 10 7 M , with 90% 
of the sample in the range 10 7 — 10 9 M Q . 

Figure 16 (left) shows the calculated dust 
masses for the SSGSS sample versus SFR e , with 
colour coding to indicate AGN type by BPT des- 
ignation. Also plotted are the results of |da Cunha| 
eTHl pHOt (D10), who fit the GALEX-SDSS- 



2MASS-/iL4S' photometry of SDSS and SINGS 
galaxies with SEDs derived from population syn- 



thesis and dust attenuation models (da Cunha, 
|Charlot fc Elbaz] |2008| ) to calculate both dust 
mass and SFR. SSGSS dust masses are very sim- 



(2007 


) models for 


ies ( 


Draine et al. 



comparibility of the two models. 

Mdust is highly correlated with SFR e , with an 
almost linear relationship that is similar to that 
found by D10. The best-fit regression line is: 



M dust (M ) = 
(2.0 ± 0.3) x 10 7 SFi^MoZjr- 1 / - 95 * - 11 ) (2) 

We note that D10 find a tighter relation be- 
tween these two parameters for SDSS galaxies. 
This may indicate a more accurate SFR estimation 
by D10 over the optical emission line estimate, due 
to their use of the full FUV-to-FIR SED. 



Figure 16 (middle) shows dust-to-stellar mass 
fraction versus SSFR for SSGSS, again including 
the results of D10. There is a significant posi- 
tive correlation between dust fraction and SSFR 
(t = 0.24; correlation significance: p — 0.001), 
although here there is significantly more scatter 
than observed by D10 in their sample. The best- 
fit regression line is: 



(2.5±1.0)xl0 5 SSFR(yr- 1 ) ( - a - 8±a - 2 '> 

(3) 

Figure 16 (right) shows the ratio of dust mass 
to SFR, which, given the tight correlation between 
SFR and gas mass, may be considered a proxy for 
the dust-to-gas mass ratio, and so traces the en- 
richment of the interstellar medium. The SSGSS 
results follow the general trend of the D10 sample, 
(r = —0.175; correlation significance: p — 0.015), 
however the careful sample selection of SSGSS 
means that it contains fewer of the rare, extreme 
SF galaxies, and the UV selection eliminates the 
oldest red, dead galaxies; the result is a more re- 
stricted range of SSFR that makes it difficult to 
verify the slope of this relation. 

These strong correlations between galaxys' dust 
content and star formation rate — and especially 
between their absolute values — supports the idea 
that SFR provides a useful proxy for dust mass, 
as found by D10. The careful sample selection of 
SSGSS adds weight to this result. 
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Fig. 16. — Left: dust masses for SSGSS galaxies from models of Draine & Li (20071 versus SFR e (Brinchmann 



et al. 2004) (coloured circles, coded by BPT designation; see Fig. 13). Also plotted are the results of the 
dust model of da Cunha et al. (2010) for SDSS galaxies (black points) and SINGS galaxies (pink crosses). 
The best log-log regression lines are shown for SSGSS (dashed) and SDSS (dotted). Middle: dust-to-stellar 
mass fraction versus SSFR, and right: dust mass-to-SFR ratio. 
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Fig. 17. — Upper Panels: dust-to-stellar mass fraction versus the optical emission line ratios [OIII]/H/3 
(left) and [Nil] /Ha (middle). These line ratios define BPT AGN type, which is designated by the same 



colour coding as Fig. 13 Also shown is the distribution of dust fractions by AGN type (right). AGN have 



significantly lower dust fractions than SF galaxies. SF galaxies with [OIII]/H/3 > 0.5 (dotted line) are young 
galaxies with starburst-source hard radiation fields; these galaxies have significantly higher dust fractions 
than AGN with similar [OIII]/H/3. Lower Panels: same as above, except with the ratio of dust mass to 
SFR P . 
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As noted by D10, these relations can be inter- 
preted in evolutionary terms. The ISM is enriched 
by dust with ongoing star formation, and at the 
same time gas mass decreases, resulting in a drop 
in SFR. With this drop, the production of dust 
is quickly superseded by destruction of grains in 
the ISM, resulting in the positive correlation be- 



tween dust fraction and SSFR (Fig. 16 middle) 



At the same time, this reduction in dust content 
with galaxy age appears to be outpaced by the 
reduction in gas content and the associated drop 
in SSFR, resulting in the anticorrelation between 
dust-to-gas mass ratio vs. SSFR (Fig. 16 right). 

In this scenario, the balance between ISM en- 
richment and the destruction of dust by shocks 
and by the ambient radiation field defines the di- 
rection and slope of these relations. The BPT 
emission line ratios, [OIII]/H/3 and [Nil] /Ha, are 
sensitive to the hardness of the ambient radiation 
field and the level of ISM enrichment respectively, 
and so may shed light on the nature of this bal- 
ance. Additionally, as seen in Figure 16, galax- 



ies with BPT-designatcd AGN components have 
lower dust fractions than SF galaxies. Although 
this may be purely due to the increased incidence 
(and observability) of AGN in older galaxies, the 
difference in dust fraction with SSFR seems to be 
divided surprisingly cleanly down AGN lines. It 
is possible that AGN shocks and hard UV pho- 
tons are important in destroying small dust grains, 
given their potential affect on PAH grain size dis- 
tributions ( |Smith et al.]2007| OD09; T10) 



Figure 17 (upper panels) shows dust-to-stellar 
mass fraction versus [OIII]/H/3 and [Nil] /Ha, as 
well as the distribution of dust fractions by BPT 
designation. Both ratios show a scattered but sig- 
nificant inverse correlation, and in both cases there 
is no significant trend taking either SF galaxies or 
AGN separately. 

[OIIIJ/H/3 provides a good proxy for the hard- 
ness of the UV radiation field, and so it is tempting 
to link dust grain destruction with an increase in 
hard UV photons, although it is not possible to 
extricate this effect from the drop in dust produc- 
tion with diminishing star formation. Notably, SF 
galaxies with hard UV fields ([OIIIJ/H/3 > 0.5) all 
have have significantly higher dust fractions than 
AGN with similar [OIII]/H/3. These same galaxies 
have among the highest SSFRs, implying that ei- 
ther the radiation field has little to do with grain 



destruction, or that dust production due to heavy 
star formation in these galaxies more than com- 
pensates for such depletion. 

[Nil] /Ha is sensitive to gas-phase metallicity, 
and so may be expected to have a positive cor- 
relation with dust-to-stellar mass fraction. The 
observed negative correlation supports the idea 
that the drop in dust content of these galaxies 
(Md us t/M*) outpaces the increase in the gas-phase 
metallicity (and corresponding dust-to-gas mass 
fraction) as these galaxies evolve. 



Figure 17 (lower panels) shows the ratio of dust 
mass to SFR vs. BPT line ratios. There is no cor- 
relation between this and any of the line ratios, 
and no trend with AGN status. This is expected; 
while dust-to-stellar mass fraction decreases with 
stellar population age, gas content and the associ- 
ated SFR probably decreases more quickly, elimi- 
nating any trend. 

Dust masses for the entire SSGSS sample are 



available as part of the data release (see Sect. 8.2) 



8. Data Products and Value-Added Cata- 
logs 

8.1. Spitzer Data Archive 

SSGSS is a Spitzer Legacy Program, and all 
spectroscopy is available for download from the 
Spitzer Data Archive. The archive is accessed via 
Spitzer's Leopard software. The web interface for 
the archive and link to the Leopard download page 



can be found at: http://archive.spitzer.caltech.edu/ 

The data release includes low-resolution spec- 
tra of the 100 galaxies in the full SSGSS sample 
and high-resolution spectra of the 33 galaxies in 
the bright SSGSS sample. These are reduced as 
described in Section [5j We include spectra with 
slit positions both combined and uncombined, and 
with orders and (for low-res) SL and LL modules 
both stitched and unstitched. 

8.2. Complementary Data Products 

In addition to the spectroscopic data, a range 
of complementary data products are available via 
the SSGSS website: 



http://astro.columbia.edu/ssgss 



UV through IR imaging data: including FUV and 



27 



NUV from GALEX; ugriz from SDSS; IRAC 3.6, 
4.5, 5.8, and 7.8 fim bands; MIPS 24, 70, and 
160 /im bands; and 16 /im IRS Peakup from SS- 
GSS, with corresponding uncertainty and/or back- 
ground images. 

Matched catalogs of SDSS properties: these are 
from the SDSS studies at MPA/JHU grourQ 
( |Brinchmann et al. 2004), and are drawn from 
data release 4 of their catalogs. These include 
measured and derived parameters such as magni- 
tudes, line strengths, Petrosian radii, SFRs, and 
index strengths, as well as magnitudes K-corrected 



using the method of Blanton et al. (2006), stel 



lar masses, mass-to-light ratios, dust attenuations, 
etc. from Kauffmann et al. (2003a), and gas-phase 



metallicities from Tremonti et al. (2004) 



PAH strengths: integrated intensities for the most 
prominent PAH features and complexes at 6.2 /im, 
7.7 /im, 8.6 /im, 11.3 /im, and 17 /im, as well as 
total PAH intensity, measured by the PAHFIT 



(Smith et al. 2007) decompositions of SSGSS spec- 



tra as described in OD09. 

Emission line strengths: line strengths from both 
the low-res and hi-res spectra from PAHFIT de- 
compositions. 



Derived properties: dust masses (see Sect. 7.4), 



silicate strength (see Sect. 7.2), and silicate opti- 



cal depth from PAHFIT decompositions. 



9. Summary 

We have presented SSGSS, a new Spitzer spec- 
troscopic survey of galaxies, colour-selected to be 
representative of the normal, star-forming galaxy 
population. SSGSS targets 101 galaxies in the 
Lockman Hole, all with extensive ancilliary imag- 
ing and spectroscopy from the FUV to the FIR. 
All galaxies were observed with IRS low resolution 
modules from 5 — 40 /^m and a subsample of the 
33 brightest were observed with at high resolution 
from 10 — 20 /im. 

For the low resolution spectra, our target S/N 
of 5 per resolution element was exceeded except 



in some SL2 spectra and in regions of high silicate 
absorption around 10/im. At high resolution, S/N 
> 10 was achieved for all but one [Ne II] lines and 
S/N > 10 for -90% of [Nc III] lines. Of the 101 
observed galaxies, 100 resulted in useful spectra, 
which are available from the Spitzer archive. 

Comparison of SSGSS galaxies to starbursts 
shows that the starburst spectra differ primarily 
in the shape of the warm dust continuum, consis- 
tent with having elevated emission at small-grain 
dust temperatures of < 120 K. Silicate absorption 
levels of SSGSS galaxies are uniformly lower than 
those often observed in LIRGs and ULIRGs, and 
are similar to those observed in starbursts. 

Using PAH strengths measured with the PAH- 



FIT code of Smith et al. (2007), we have inves- 



tigated the link between total PAH emission and 
both IR luminosity and star formation indicators. 
PAH luminosity correlates with star formation in- 
dicators -D„(4000) and Ha EW very tightly, except 
for the oldest ~20% of the sample. PAH-to-total 
IR luminosity ratio changes slightly between SF 
galaxies and those with AGN components, drop- 
ping from —0.07 to —0.05. This approximately 
reflects the change between actively star forming 
galaxies and more evolved galaxies. 

Dust masses calculated from the models of 



Draine & Li (2007) reveal a close relationship 



between dust mass and SFR, and are consistent 
with the evolutionary scenario described by D10 
in which dust that is produced quickly in a rapidly 
star-forming phase is destroyed more quickly than 
it is replenished as star formation drops. 

SSGSS offer great scope for further study. In 
OD09, the sample has been used to conduct a de- 
tailed analysis of the physical drivers of PAH spec- 



tral shapes, and Treyer et al. (2010) investigates 



1 see www. mpa-gar ching .mpg.de/SDSS/ 



the use of MIR diagnostics for SF and AGN activ- 
ity. Future work will include the development of 
template PAH spectra calibrated to a wide range 
of observables, analysis of the role of the size, 
shape, and composition (via metallicity) of small 
dust grains in determining the continuum shape. 
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Appendix A: Multi- Wavelength Thumbnails 
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Fig. 18. — Multi-wavelength images of the 101 galaxies of the SSGSS sample. Each horizontal strip shows 
a given galaxy in, respectively: GALEX FUV and NUV, SDSS multi-colour, r, H, IRAC 3.6 /im and 8 /im, 
16 /mi from SSGSS peakup images, and MIPS 25 /im and 70 /im. All boxes are 60" on a side except for the 
70 /im band, which is 120". 
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Appendix B: Reduced IRS Spectra 
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Fig. 19. — Short-low and long-low IRS spectra, reduced and stitched as described in Section [O] The yellow 
shaded region shows the ±l-sigma errors. The wavelength range of each IRS module is shown at the base 
of each panel, from left to right: SL2 (blue), SL1 (red), LL2 (blue), LL1 (red). Also given are the SSGSS 



catalog number, BPT type (Baldwin et al. 
diagnostic. 



19811, redshift, stellar mass, NUV-R color, and £>„(4000) age 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 



50 



i 74 Composite Source 
1 fc" z-O.ll NUV-R-a.flO 
H,= 11. 33 D„(4[MH)=1-J7 




! 78 Star Forming Galaxy 
z-fl.17 NUY-R-2.74 
K,= 10.3S D,(4n00)=L.Z2 



i 78 Star Forming Galaxy 
1 ^ z-G_£l74 NUV-R-l.a3 
«j 9.295 11,(4000)= 1. 13 



i 77 Star Forming Galaxy 
z-fl_fl41 HW-tt-tM 
K.= 10.4a D„(4[K»}=1.ZB 




55CSS ?B SLar Forming Galaxy 
z = 0.D74 nu*-li = S.ol 
U a =B.B32 D,(4DOa)=l.Z7 



SSGSS 79 Star Forming Galaxy 

z-a.au NiTf-R-ua 

H„= 10.47 D,[4[]0U}=I.ZS 



SSGSS A) Star Forming Galaxy 
1 |" z = D.D79 HUV-R=S.!!0 

io 22 n,(4Doa)=i.3y 



SSGSS Si Star Forming Galaxy 

z=o.07S nuy-r=2.s7 

«„= 10.74 D,[4[]0Q)=l.ZI 




30 8 10 
Rest Wavelength (^m) 



Figure 19 (continued) 
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Figure 19 (continued) 
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Figure 19 (continued) 
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Fig. 20. — Short-high IRS spectra, reduced and stitched as described in Section 5.2 The yellow shaded 



region shows the ±l-sigma errors. Also given are the SSGSS catalog number, BPT type (Baldwin et al 



1981), redshift, stellar mass, NUV-R color, and Z?„(4000) age diagnostic. 
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Figure 20 (continued) 
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Figure 20 (continued) 
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Figure 20 (continued) 
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Figure 20 (continued) 
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